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We show the equivalence of two approaches to the design of entanglement distillation protocols. 
The first approach is based on local unitary operations that yield permutations of tensor products 
of Bell states. The second approach is based on stabilizer codes. 
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I. INTRODUCTION 

We describe a link between two approaches to the de- 
sign of entanglement distillation protocols. The first ap- 
proach, which we call permutation based, has it origins in 
Refs. [li!3], where protocols are presented based on local 
operations that permute the 16 possible tensor products 
of 2 Bell states. This approach was generalized in Ref. .^ 
through a description in binary arithmetic of all locally 
realizable permutations of all tensor products of n Bell 
states. Several other protocols based on similar ideas as 
in Refs. JmH have been proposed. We only mention 
Refs. 0, H H H- The second approach, which we call 
code based, derives distillation protocols from stabilizer 
codes as in Refs. @,l3- 

In the permutation based approach the protocol con- 
sists of unitary operations applied by Alice and Bob to 
their side of a number of qubit pairs (in some of the ref- 
erences above also higher dimensional systems are con- 
sidered, but in this paper we will restrict ourselves to 
qubits) , followed by a number of measurements on some 
of the pairs (single qubit measurements) . The measure- 
ment outcomes are compared to decide if the protocol 
has succeeded. If two-way communication is possible, the 
protocol can be followed by new stages, possibly affecting 
more qubit pairs. 

In the description of code based protocols a number of 
commuting joint measurements is performed by Alice and 
Bob (corresponding to measuring the syndrome of a sta- 
bilizer code). The measurement outcomes are compared 
in the same way as for the permutation based protocols. 
To end up with Bell states a final operation is needed at 
both sides that corresponds to the decoding operation of 
the stabihzer code. 

We show that both approaches are in correspondence 
in the following way. If the measurement of the code 
based protocol is performed by first decoding, then do- 
ing qubit measurements and then encoding again, the 
protocol becomes equivalent to a permutation based pro- 
tocol, as the encoding operation is canceled by the final 
decoding operation of the code based protocol. The de- 
coding operator of a stabilizer code applied at both sides 



yields a permutation of tensor products of Bell states. 
Conversely, every permutation based protocol can be in- 
terpreted as a code based protocol. 

In this paper we will focus on two schemes presented 
in Refs. 0jl3 that represent the two approaches. In sec- 
tion ^ we give an overview of some preliminary defini- 
tions and theorems. In section UTTl we describe the main 
results of the permutation based and the code based pro- 
tocols, omitting the more detailed elaborations, as stated 
in Refs. y,|^. Both of them are slightly moderated in 
order to clearly show the equivalence, which is done in 
section Hvl 



II. PRELIMINARIES 

A. Binary vector representation of products of Bell 

states 

Bell states can be represented by assigning two-bit vec- 
tors to the Bell states as follows 
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A tensor product of n Bell states can then be described 
by a 2n-bit vector, e.g. |Boiooii) = |-Boo)|-Bii)|-Boi) = 
|$+)|$-)|ijf+). Note that the first half of the subscript 
refers to the phase of the Bell states ("-I-" or "— "), and 
the second to |$) or |^). 

An interesting feature is the correspondence between 
Bell states and Pauli matrices 
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A tensor product of n Bell states is then de- 
scribed by a kronecker product of Pauli matrices, e.g. 

|$+)|*-)|*+) -^ -js'^o ® CTy ^ ax == TT^croO «> crii 0(701- 
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In the same way we can use longer vector subscripts to 
denote such kronecker products, e.g. CToiooh = "'oo '^ 
an ctqi. All the 4.4" tensor products of Pauh matrices 
multiplied with a global phase £ {±1, ±i} form the Pauli 
group Vn- 

The Pauli group and the group of binary vectors Zj", 
are related to each other in the following way 



^a^b ^ ^a+b, 



(3) 



The maximum value pa of the px in © can be used as 
an entanglement measure. Note that one can always 
transform po into this maximum value by applying at 
either Alice's or Bob's side the corresponding Pauli oper- 
ator (Ta that rotates \Ba){Ba\ into \Bo){Bo\. We will call 
Po = {Bo\p\Bo) the fidelity F. 



Stabilizer codes 



where "'^" denotes equality up to a global phase [lOl| . 
Such a phase is irrelevant when these matrices represent 
pure state vectors. All addition of binary objects is done 
modulo 2. 

Note that any two elements of the Pauli group either 
commute or anticommute. It can be easily verified that 
two Pauli matrices CTq and cTf, commute if the symplectic 
inner product a^ Pb is equal to zero, or 



o-ao-fe = (-1)" abaa, where P = 
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B. Local permutations of products of Bell states 

In general a pure bipartite state lip) of 2n qubits can be 
represented by a 2" x 2"-matrix ^ (e.g. as in (0) for Bell 
states). Local unitary operations jV') — + (Ua ® C/_b)|i/'), 
in which Alice acts on her n qubits (jointly) with an 
operation Ua and Bob on his n qubits with an operation 
Ub , are then represented by 



^ 
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Interesting local unitary operations for the protocols un- 
der consideration are those that result in a permutation 
of the 4" tensor products of n Bell states (up to a global 
phase). In Ref. 3] it is proven that if local unitary op- 
erations result in a permutation of the products of Bell 
states, this permutation can be represented in the binary 
vector representation as an affine symplectic operation 
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with Aell'^''^'\b£'Ef' 
and 



A^PA. 
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and that any such permutation (j) can be realized by local 
unitary operations 



^x — > Ua^xUb = (TAxA 



(7) 



For an efficient way of doing this (by means of 0(n^ one 
and two-qubit operations), we refer to Refs. IS, [lfl|. A 
matrix satifsying © is called P-orthogonal or symplectic. 
The examined protocols result in a new mixture of ten- 
sor products of Bell states 
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A stabilizer S" is a commutative subgroup of the Pauli 
group Vn which does not contain —1 or ±j. The stabilizer 
code Tis C 7^21 associated with S is the joint eigenspace 
of all elements of S, or 
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The group S is called the stabilizer of the code, since it 
preserves all of the states in the code. If the stabilizer 
has n^m independent generators Mi , the code space has 
dimension 2™ [Ij. From the definition of S it follows 
that the eigenvalues of each Mi are -1-1 and —1, with the 
same multiplicity. 

Suppose an error E E Vn has afflicted a state |V) G Hs- 
If E and Mi commute, then 



M,E\ij) ^ EM,\^) ^ E\ib), 

so the error preserves the value -t-1 of Af, 
anticommute, then 



(10) 



If E and Mi 



M,E\ib) = ~EM,\^) = ~E\ij), 



(11) 



so that the error flips the value of Mi , and the error can 
be detected by measuring Mi. Repeating this procedure 
for every generator Mi of S, we may write 



M,E={-IY^EM,. 



(12) 



The Si, i — I, . . . ,n — m constitute a syndrome for the 
error E as (—1)"' will be the result of measuring Mi if the 
error E has occurred. For recovery, a Pauli operator R 
is applied that has the same commutation relations with 
the Mi as E. Recovery is successful if RE € S, in which 
case error -I- recovery has a trivial effect on states E Hs- 
Measuring the syndrome s comes down to projecting 
onto the joint eigenspace of the Mi with eigenvalues Si. 
Note that this eigenspace is the code space associated 
with the stabilizer with generators (— l)**Mi. We will 
call this code space Hs,s- 



III. TWO SCHEMES FOR CREATING 

ENTANGLEMENT DISTILLATION PROTOCOLS 

A. Permutation based protocols 

A slightly generalized variant of the distillation proto- 
cols presented in Ref. 3 can be summarized as follows. 



1 . Start from a mixture of 4" tensor products of Bell 
states. Typically, this is the tensor product of n 
identical independent Bell diagonal states. 

2. Apply a permutation of these 4" products of Bell 
states with local unitary transformations as de- 
scribed in the preceding section. 

3. Check whether the last n — m qubit pairs are |$) 
or |^)-states. This can be accomplished locally by 
measuring both qubits of each pair in the 1 0) , 1 1 ) ba- 
sis, and checking whether both measurements yield 
the same or the opposite result. 

4. Perform m single-qubit Pauli operations to Bob's 
remaining qubits as described further. This comes 
down to rotating the resulting state ((SJ where Pa is 
maximal to a state where Pq is maximal. 

5. If the resulting state of the m remaining pairs sat- 
isfies a certain criterion (e.g. the fidelity exceeds 
a certain proposed value), keep them, otherwise, 
discard them. The result is a new mixture of 4™ 
products of Bell states. 

Using the same techniques as in Ref . _3] , one obtains the 
following. If Alice and Bob apply the above protocol, 
starting from an initial state 



Y. P.\B,){B,\ 



(13) 



x& 



the resulting state of the m remaining pairs after the 
measurement is 



E 



yGZ 



E 



x£S+PA^Py 



r.Px 



Z^x^S^+PATPO Px 



\By){By\ (14) 



where S is the subspace spanned by the rows of AP with 
indices n + m + \, . . . ,2n and iS^ is the subspace of all 
the binary vectors that have a symplectic inner product 
equal to zero with the elements of S. Note that S C S^ . 
y € Z|" is constructed from y € Zj™ as follows 



y = yiy2 . . . y™ po . . . qyn+iyn+2 ■ ■ • yri+m tit 
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(15) 
where t is the outcome of the measurements of the last 
n — m qubit pairs, i.e. zeros where the measurements are 
the same, ones where they are opposite. 

In step^l Bob applies to his remaining m qubits the Ua 
for which the coefficient of \Ba){Ba\ in H14|) is maximal, 
or 



B. Code based protocols 

The protocols described in Ref. Q are derived from 
a stabilizer S with generators Ug^ for i — 1, . . . ,n ~ m. 
In Ref. jS|, the protocols are more generally defined for 
higher-dimensional systems (qudits). We only consider 
qubits here. The protocol goes as follows: 

1. Alice measures the observable a* for each i, where 
"*" stands for elementwise complex conjugation. 
Let (—1)°' be the results of the measurements: 
they are the eigenvalues of the stabilizer code space 
'Hs' ,a containing the state after measurement. 

2. Bob measures the observable Cg. for each i. Let 
(—1)''' be the eigenvalues of the stabihzer code 
space Ti-s.b containing the state after measurement. 

3. Bob performs the error correcting process treating 
the string s = (6i -I- oi, . . . , hn-m + an-m) (mod 2) 
as a syndrome as described below. 

4. Alice applies to her share of the qubit pairs the 
inverse of the encoding operator of Ti.s',a- Bob ap- 
plies the inverse of the encoding operator of Ti-s.b 
to his share. 

5. They discard the last n — m qubit pairs. 

6. If the fidelity of the resulting state (which depends 
on the difference of the measurement outcomes 
(6i -I- ai, . . . , bn-m + fln-m)) of the m remaining 
pairs is large enough, they are kept, otherwise, dis- 
carded. 

In stepini Bob has to apply a recovery operator R treating 
the string s as a syndrome. R is defined as follows. Let C 
be the hnear subspace of Z|" generated by 51, ... , g„_m, 
and C^ the orthogonal space of C with respect to the 
symplectic inner product. Choose v G Z2" such that a^ 
has commutation relations with the CTg. corresponding to 
s, or v^Pgi = Si. Now define u as the vector having 
maximum 
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in the coset C^ + v. Maximum fidelity (after step EJ is 
achieved if Bob applies i? = cr„ in step O and is equal to 
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a = argmax 
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Note that a solution is not unique: any element in the 
coset a -I- iS will do. The resulting fidelity is then 



F = 2" 



l^xeS+PA^PgPx 
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IV. EQUIVALENCE OF THE TWO 
APPROACHES 



We will first show that measuring the observables in 
the code based protocols can be carried out by local 
unitary operations that yield a permutation of the ten- 
sor products of the Bell states, followed by single qubit 



measurements and the inverse of the local operations. 
Second, we show more specifically that the protocols of 
Refs. and Q have the same results, i.e. the formulas 
for the resulting fidelity in Refs. Q and [g are equivalent. 



A. Observable measurements in the code based 
protocols 

Let us first consider Bob's actions. A way of measuring 
the observables Cg. is by first applying U\ which is the 
inverse of the unitary coding operator U , then measur- 
ing the last n — m qubits in the |0), |1) basis, and then 
applying U. For every b G Z^""^, U transforms states 



|00..J)6i62...6„_™),. 
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^6162... 6„-m) (20) 



to a basis of the stabilizer code space Hs,b- U satisfies 
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has 1 in position m + i and zeros 

^ J0im+i-l) (g, (J (g, J0in-m-i)^ ^ 



where e. 

elsewhere, or ae 

way of efficiently implementing such coding operator U 

is explained in Ref. [T^l- The U in (|21|l and the inverse 

of the coding operator W in step 01 of the protocol cancel 

each other, so Bob only applies U^ and then measures 

his last n — m qubits. 

In the same way, we can see that, in order to measure 
the observables cr* , Alice applies U'^ and then measures 



her last n — m qubits. 
(|2I1 it follows that 
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(22) 



If we interpret the a in this equation as representing ten- 
sor products of Bell states, we see with © and Q that 
the local unitary transformations Ua and Ub yield in the 
binary picture a permutation (f>, with 
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Vi. 



(23) 



Note that a*„ 



(T„. and these matrices represent the 



same pure states as a global phase is irrelevant. 



B. Equivalence of the t^wo schemes 

To point out that formulas IjlTI) and (|19|l are equiva- 
lent, we have to show that the subspaces S and C are the 
same and that the summations are done over the same 
cosets of 5 = C and S^ = C^ . 

We define a symplectic matrix B with Bcm+i — 
gi, which means that colums m + 1 through n are 
gi, . . . ,gn-m- A matrix with these columns can al- 
ways be completed to a symplectic matrix since S is 
a commutative group, and therefore we already have 



that gfPgj — 0,Vi,j. The inverse of this matrix is 
B^^ — PB P, which is also symplectic. B~^ defines the 
permutation realized by Ua and Ub- Note that fixing B 
by adding columns to B comes down to fixing bases of 
the stabilizer codes 'Hs*,a and Ti-sfi- The last n — m rows 
of B^^P = PB^ are 51, ... , gn~m, which means that S 
and C in H17() and (|19|l are the same. 

We have C^ + PP-^PO ^C^ +vii PP-^PO = BO 
has symplectic inner products with gi, . . . , gn-m equal 
to si, . . . , Sn-m- We know from the definition of s and t 
that s — t. We then have 



1=1 



(24) 



Z^" 



with Bk the fc-th column of B. We know from the sym- 
plecticity of B that Bn+m+i has symplectic inner prod- 
uct equal to zero with all the other columns of B except 
i3m_l_i = gi, from which follows that {BO)^Pgi = Si 

In the same way it can be verified that Vy G 
By e C^ + BO, OT C + Ba C C-^ + BO so the cosets 
C + Ba and C + u are the same since they are both 
cosets of C in C^ having maximum total probability of 
their elements. 



V. CONCLUSION 

We have compared two approaches to the design of 
entanglement distillation protocols. Although both ap- 
proaches are based on different concepts, they turn out 
to be to a large extent equivalent. 
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